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We have performed molecular dynamic simulations of a Martensitic 
bcc-^hcp transformation in a homogeneous system. The system evolves into 
three Martensitic variants, sharing a common nearest neighbor vector along a 
bcc (111) direction, plus an fee region. Nucleation occurs locally, followed by 
subsequent growth. We monitor the time-dependent scattering S'(q, t) during 
the transformation, and find anomalous, Brillouin zone-dependent scattering 
similar to that observed experimentally in a number of systems above the 
transformation temperature. This scattering is shown to be related to the 
elastic strain associated with the transformation, and is not directly related 
to the phonon response. 
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I. INTRODUCTION. 

The topic of precursor phenomena in martensitic transformations has been the subject of 
numerous studies over the past two decades. Martensitic systems undergo first-order tran- 
sitions which are characterized by a lattice strain (sometimes accompanied by a "shuffling" 
of the atoms associated with a particular phonon mode). 00 Experimentally, it has been 
observed from various diffraction measurements that streaking of the Bragg spots occur in 
the parent phase close to the transformation, indicating an incipient instability towards the 
formation of the martensitic phase.H^ There is also evidence that interesting anomalous dy- 
namical behavior occurs concurrently: the experimentally observed phonon peaks become 
very broad with enhanced scattering occurring at low frequencies.oQ More puzzling, the in- 
tensities and shapes of these peaks exhibit unexpected variations between different Brillouin 
zones. While most of the experimental data were obtained from alloy systems where the 
effects of stoichiometry variations may be of importance, similar behavior have also been 
observed in pure metals.iHii 

In this paper, we explore the relationship between anomalous phonon behavior and the 
phase transformation process through molecular dynamics simulations of Zr as a prototypical 
system, exhibiting a bcc to hep transformation. By monitoring the behavior of the dynamic 
structure factor of the system, we examine the microscopic origin of the anomalous phonon 
behavior in these systems. We also study the dynamics of microstructure formation during 
the transformation. To clarify the problem, we will concentrate on homogeneous systems 
without defects. The effects of inhomogeneity on the dynamical behavior of the transition 
will be left to a later study. 
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Scattering lineshapes in first-order structural phase transformations liave been previously 
studied using various analytic treatmentJiilli and simulationJlUlZl. Such studies, for both 
first-order and second-order transformations, have focused on the one-phonon contribution to 
the scattering, which does not exhibit a Brillouin zone dependence, even in a full anharmonic 
treatment. Furthermore, most of these studieJii^ have focussed on models that do not 
exhibit a true elastic transformation: in the high temperature phase, the linearized phonon 
frequencies behave as uj"^ ^ uJq + c^k"^ as — > 0, rather than uj — > ck. Such models are more 
appropriate for transformations such as hcc uj, characterized by a "freezing in" of a phonon 
mode with a finite frequency (and wavevector), rather than the development of a static 
amplitude of an elastic degree of freedom that occurs in many Martensitic transformations 
(including the hcc — hep transformation). 

Another issue with previous work examining the dynamic response of these systems is 
that most work examines regions of the product (low temperature) phase embedded in the 
high temperature matrixlli'@lli. This product phase is hypothesized to be due to either 
defects0'ti, or due to long-lived "heterophase" fluctuations from the parent phase. Such 
studies are not surprising, as long-lived excitations are necessary to produce "central peaks" 
in the scatteringEl. However, in first-order transformations, the evidence appears to indicate 
that there is no clear central peak separate from the phonon peaks (and no long-lived 
fluctuations into the product phase), especially in systems with little disorder. 

In a previous paper ,0 which we shall refer to as paper I, we demonstrated that in a homo- 
geneous system, fluctuations in the high-temperature hcc phase can produce "anomalous" 
scattering, with significantly different scattering in equivalent Brillouin zones. Neutron 
scattering reveals different line shapes between different Brillouin zones.ilfl'i Thermal dif- 
fuse scattering (energy integrated scattering), observed in both X-ray studies! and neutron 
diffraction experiments^, is not featureless, as might be expected from usual arguments 
concerning multiple phonon scattering,!^ but instead shows characteristics that imply a con- 
nection with the incipient Martensitic transition. That such characteristic scattering can 
occur in the absence of static or long-lived fluctuations into the product phase is somewhat 
surprising, yet is borne out by experimental work. Hi In paper I, we argued that these are 
due to fluctuations toward, but not into, the low temperature phase. 

More recent molecular dynamics simulationf3cJ have examined the formation of mi- 
crostructure during a homogeneous hcc-^hcp transformation, utilizing a simple EAM model 
for zirconium, as well as the effect of stress on the microstructure. These simulations demon- 
strate a number of key points concerning the formation of microstructure. Of particular 
importance relevant to the work presented here, they have demonstrated that fee regions 
form at the junction of three hep variants sharing a common (1120) direction, rather than 
forming a triple point between twin boundaries, and that the twin regions contain basal 
plane stacking faults that maintain the 60° orientation between hep variants. The transfor- 
mation is interpreted to occur via the Nishiyama-Wassermann mechanismll (described in 
Section III below), rather than the Burgers mechanism.il 

Here, we present similar simulations, with the goal of examining the dynamics during the 
hee-^hep transition, with a focus on the aspects relevant to the neutron and x-ray scattering 
experiments. Again, we use molecular dynamics as a tool, to examine the microscopic pro- 
cesses associated with the dynamics. We examine the dynamics in real space, to examine the 
growth process and the final microstructure. The crystallography of the resulting twinned 
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structure corresponds to experimental results, and those of PinsookEra. With respect to 
the microstructural development, the most significant difference between the work presented 
here and that of Ref. is the choice of simulation cell orientation: unlike their simulations, 
we do not select out one of the bcc directions to be alone one of our simulation that 
favors hep variants with (1120) along this direction. 

We also examine the time-dependent structure factor S{Q,t), which is the equivalent of 
time-dependent x-ray scattering during the transition. By doing this, we directly demon- 
strate that the asymmetry in the scattering between different Brillouin zones is strongest 
during the transition, and is essentially due to the formation and motion of new Bragg peaks 
during the transition. 



II. DETAILS OF SIMULATIONS. 

We performed the simulations using the same EAM potential as in the previous paper ,0 
and described in detail in Ref. ^ This potential has a stable bcc phase, with phonons 
in good agreement with experimental results. The difference between the cohesive energies 
between the fee and the hep structures is small, compared with ab initio results, as is common 
for these types of potentials. The EAM potentials for Ti and Zr are prone to this problem, 
due to the lack of directional bonding in these descriptions of these materials. Another 
important caveat is that this potential does not accurately reproduce the w-phase properties 
of Zr or Ti: within this potential, the w-phase structure is mechanically unstable. As a result, 
dynamics associated with large-amplitude fluctuations towards this structure (corresponding 
to the bcc q=|[lll] longitudinal phonon) are not accurately reproduced. In Zr and Ti (and 
especially in certain alloys, including ZrNb), the large amplitude fluctuations of this will be 
anharmonic and possibly strongly non-linear, associated with the corresponding bcc — >■ u 
transformation. This latter transition, while not studied in this work, is an important 
topic. The work presented here is essentially a model system, examining the dynamics of 
the bcc— >■ hep transformation in the simple case where no competing transformation occurs. 

The simulations were performed using two different cubic system sizes: 18 x 18 x 18 
conventional bcc unit cells (11664 atoms), and 36 x 36 x 36 unit cells (93312 atoms). Except 
as noted, all results presented will be those of the larger system. No quahtative changes 
occurred as a result of the change in the system size. Periodic boundary conditions were 
used in all cases. In order to minimally interfere with the dynamics during the transition, the 
simulations were performed under constant energy conditions, at an energy corresponding to 
the average energy for the bcc phase at 1500 K. During the bee— >■ hep transition, the system 
lowered its potential energy, as can be seen in flg. |^. With constant energy simulations, 
this produces a slight increase in temperature. For both the 18 x 18 x 18 system and the 
36 X 36 X 36 system, this increase was about 40 K. This is a minor change in temperature, and 
we do not believe that this affects the dynamics of the transition signiflcantly. We have also 
used constant volume simulations, with a bee lattice constant of 3.68 A as in paper I. The 
transformed structure does not develop any signiflcant macroscopic stress. The simulations 
were performed for 40,000 time steps, sufficient to observe the transition and subsequent 
equilibration. This will be demonstrated in the following section. 
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III. DEVELOPMENT OF MICROSTRUCTURE DURING THE TRANSITION. 



We begin with a review of the crystallography of the bcc-^ hep transition. This is most 
commonly stated in terms of the Burgers relationship,ii 

(110)5ec||(0001),,, [lll],,,||[1120];,,p. (1) 

The transformation occurs by "shuffling" of alternate (110) planes in the bcc lattice, simul- 
taneous with a shearing of the (112) planes along the [111] direction. (In this paper, we will 
commonly use three indices for directions and planes in the bcc lattice, and four for directions 
and planes in the hep lattice. We will indicate the lattice structure explicitly where we wish 
to emphasize the lattice correspondence, but otherwise will omit these labels.) The (112) 
direction then becomes a prism plane in the hep lattice. With these dynamics in mind, we 
may rewrite the Burgers relationship asil 

(112),,,||(1100);,,p; [lll]5,J|[1120];,,p. (2) 

The shearing of the lattice along {112} planes plays an important role in understanding the 
anomalous scattering, as will be shown in section 4. Note that if the [1120] /^cp direction is 
along the [lll]^^^ direction, then there are three different {110}{,cc planes that can become 
hep basal planes, namely the (lOl), (Oil) and (llO) planes. With four equivalent {lll)bcc 
directions in the bee lattice, there are 12 hep variants that may occur. In our simulated 
microstructure (described below), the system spontaneously selects out one of the (lll)6cc 
directions, and forms three hep variants from the associated { 1 10 j^cc planes. 

We now establish that there is a transition during our simulations. The most direct 
approach is to examine the potential energy vs. time during the transition. This is shown in 
fig. |l|, for both the 18 x 18 x 18 and 36 x 36 x 36 system sizes. As can be seen in the figure, 
both systems spend considerable amounts of time (10,000-20,000 time steps) fiuctuating 
about a value near -6.032 eV. This corresponds to the value obtained at 1500 K for the bee 
phase. Both systems then have a rapid drop in the potential energy. This drop continues for 
some time, then levels off near the value appropriate for the hep system. Thus, the period 
during which the potential energy changes systematically serves to approximately identify 
when the transformational dynamics occur. 

Before describing the dynamics of microstructure formation, we present the final mi- 
crostructure in the system. This is shown for the 36 x 36 x 36 system in fig. 0. Again, 
results for the smaller system are qualitatively the same. The figure shows the view down 
the [111] direction. There are three different hep variants present; the [111] direction forms 
a (1120) direction in each of these variants, in accordance with the Burgers relationship.il 
In the projection shown in fig. ^, the system appears to be very well ordered. Each of the 
variants forms a "needle-like" domain, extended indefinitely (via the periodic boundary con- 
ditions) along the [111] direction. By coloring the atoms according to their potential energy, 
we have highlighted the twin boundaries between different hep variants. These boundaries 
form along hep (jlOll) planes, in agreement with experimental observations. As in previous 
simulations,E3Ej, some of the variants have numerous basal plane stacking faults, allowing 
for the average angle between hep domains to be 60°. However, specific twin boundaries 
deviate noticeably from this value. 
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We note that in the figure, it is possible to travel from variant II to variant III without 
passing through such a twin boundary. In fact, there is an fee region formed where these 
variants intersect, as indicated in the schematic in the figure. This region is stable throughout 
the simulation, and is observed in both the 18 x 18 x 18 and 36 x 36 x 36 systems. It is 
not required to minimize the potential energy of the final configuration: a lower energy 
configuration would arise from transforming the variant III regions to variant I, and the fee 
regions to variant II. This would eliminate the total amount of twin boundary, and would 
also eliminate the unfavorable fee region. We note that the fee-hen energy difference is 
small in this potential, relative to ab-initio calculations of pure Zr,B3 suggesting that fee 
formation in Zr may be decreased or eliminated in the experimental system. However, some 
observations of fee Ti have been made in association with the bcc^hep transformation, as 
discussed in ref. 

We now examine the time evolution of the system, as it progresses from the bee to hep 
phases. This is shown for the 18 x 18 x 18 system in fig. ^. We focus on the smaller sized 
system for clarity; the results for the 36 x 36 x 36 system are quite similar. In the figure, 
we show a slice of the system perpendicular to the [110] direction. The atomic positions 
have been averaged over 2000 time steps, in order to eliminate the high frequency atomic 
motions, and reveal the slower dynamics associated with the transformation. The atoms 
have been colored according to their time-averaged potential energy. Again, this highlights 
the development of the microstructure. 

The first frame of fig. |^, shows the results for time steps 18,000-20,000. Before this, there 
is little difference between the structure and that of a perfect bcc lattice. This is consistent 
with fig. 0, which shows the potential energy of the 18 x 18 x 18 system beginning to decrease 
near step 20,000. In the first frame, the start of the transformation can be seen just left of 
the center of the system. The vertical rows of atoms show considerable bending around this 
region, indicating the large accommodation strain required for the nucleus. By the second 
frame (steps 20,000-22,000) the transforming region is quite clear. In steps 22,000-24,000, 
the transformed region is starting to develop sharp boundaries, as it extends along the [111] 
direction. This transformation progresses and the boundaries develop through step 26,000; 
by steps 26,000-28,000, the transformation is nearly complete, with only the final stages of 
boundary development remaining. After step 30,000 (not shown), the structure is essentially 
static. 

IV. SIMULATION RESULTS FOR THE STRUCTURE FACTOR. 

One of our main objectives is to understand the origin of the pre-martensitic behavior 
as seen in neutron and x-ray scattering experiments. One of the most notable anomalous 
behaviors is the difference in thermal diffuse scattering between the so-called "a;" points, with 
strong scattering located at Q = |(774), and significantly less scattering at Q = |(558). 
These correspond to equivalent points in the Brillouin zone, with equivalent magnitudes. 
Furthermore, in a simple one-phonon picture, the dependence of the scattering on the phonon 
polarizations are also identical. Therefore, from the one-phonon picture, the scattering at 
these points should be identical. These modes correspond to the |(111) longitudinal phonon, 
which is the mode that leads to the development of the cu-phase. This connection between 
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the anisotropy and possible precursors is very suggestive. However, as shown in paper I 
and elaborated on below, a significant difference between these Q vectors is related to the 
development of the hep phase, and the dynamics of the transformation to this phase. 

We begin by showing the structure factor S'(Q) calculated for a smaller system, 12 x 
12 X 12 unit cells, which does not transform. The results for scattering in the (110) plane 
are shown in fig. |. This has been calculated in the same manner as the results in paper I. 
This figure shows that within our simulations, the anomalous behavior is reproduced, with 
the scattering near Q = |(774) significantly greater than that near Q = |(558). In our 
potential, the u phase is higher in potential energy than the bcc phase, and is mechanically 
unstable. Therefore, the results we show here can not be attributed to u phase fluctuations. 
The anomalous scattering can be seen as a ridge extending from the [222] Bragg peak down 
to the I [552] N-point, where a small peak in the scattering occurs. 

We note that in these simulations of the bcc phase, there is no evidence of long-lived 
fluctuations. The pair correlation function g{r) exhibits clear peaks associated with the bcc 
phase, but none associated with the hep phase. However, some short-lived fluctuations to- 
ward the hep phase are observed, and are correlated with anomalous scattering (as described 
in paper I). This is in accordance with both experiments on the high temperature phaseiHll 
and theory.0 Thus, the anomalous scattering can not be ascribed to previously postulatedlii 
long-lived fluctuations into the hep phase. 

In order to explore the scattering behavior, we examine here the time-dependent scatter- 
ing at these "strong" and "weak" u points in the Brillouin zone. As paper I demonstrated, 
the difference in scattering between these two points can not be found within the one-phonon 
contribution, and the full scattering function must be calculated. Moreover, as there are 12 
equivalent wave vectors corresponding to each of these points, we must follow the dynamics 
of each point. 

Figure || shows S{Q,t) for Q = |(774), for the 18 x 18 x 18 system. As indicated 
previously, the transformation occurs essentially between steps 20,000 and 30,000. In the 
figure, we see that for three of the wave vectors, intense scattering occurs during the trans- 
formation. The scattering is as large as 20 times that of the other wave vectors, or of the 
same wave vectors before and after the transformation. For the other nine wave vectors, the 
scattering shows no distinctive behavior before, during, or after the transformation. Figure 
H shows S{Q,t) for Q = i(558). In this case, there is no distinctive behavior during the 
transformation. For three of the wave vectors, the scattering increases by roughly a factor 
of five once the transformation is complete. 

We can understand these results by examining the full scattering S{Q) in the (110) 
scattering plane as a function of time. Rather than perform the full calculation, we show 
the scattering calculated from the time- averaged positions, such as those shown in fig. ^ We 
show these calculations in figs. and |. Note that these figures contain essentially static 
information about the microstructure. Before the transformation is initiated (steps 14,000- 
16,000), the structure factor shows the Bragg peaks expected for the bee lattice, plus some 
weak scattering near the N-point phonons at Q = [ffl] and at Q = [§§3]. Some smearing 
of the Bragg peaks along the [112] direction is evident. In steps 16,000-18,000, the scattering 
at the N-points has intensified. 

The most dramatic change to the scattering (and the most important to understanding 
the anomalous scattering) is seen in steps 18,000-20,000. The N-point phonons are now 
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developing into static Bragg peaks, and the Bragg peaks are now shifted from their bcc lattice 
positions. The shearing along the [112] direction is evident. This shear brings scattering 
towards the Q = |(774) scattering vector, while leaving little scattering in the Q = |(558) 
direction. This is maintained through step 26,000 - essentially, until the transformation 
is nearly complete. The Bragg peaks now are clearly arranged in a hexagonal pattern. 
Thus, the anomalous large scattering near Q = |(774) is not from phonon scattering from 
I [111] phonons, but from essentially static displacements associated with N-point phonons 
combined with shears in the [112] direction. These are exactly the displacements needed 
for the formation of the hep phase. From step 26,000 on, the primary development is the 
narrowing of the peaks along the [111] direction, corresponding to the formation of long, 
needle-like domains in this direction as seen in the real space pictures in fig. ^ The Bragg 
peaks also settle into their final positions, with little remnant scattering in the Q = |(774) 
or the Q = |(558) directions. 

We also note that during the course of the transformation, a Bragg peak develops near 
Q = I [111]. This is a result of the differing hep variants (and the fee region, as well) 
developing nearest neighbor lattice vectors along the [111] direction. For real systems (i.e. 
Ti and Zr) the presence of this developing Bragg peak clouds the issue of interpreting the 
dynamic behavior of the high temperature phase at this wave vector, as anomalies can be 
associated either with the bcc —>■ hep or the bcc —>■ lo transformation. 



V. DISCUSSION. 

The main purpose of this paper has been to explore the dynamics of the bee-^hep trans- 
formation, and to relate this to the anomalous diffuse scattering observed in the bee phase. 
By examining the time-dependent scattering during the scattering, we have shown that the 
Brillouin-zone dependence of the scattering is greatest during the transformation, and is due 
to the formation and shearing of the Bragg peaks as they form the hep reciprocal lattice. The 
peaks move through the Q = |(774) region of reciprocal space, creating intense scattering 
at this wave vector, but their final position is not at this point. There is no such movement 
through the Q = |(558) point; hence the difference in scattering at these points. A one- 
phonon picture, taken from the bee lattice, would predict identical scattering at these ppinta, 
even in a fully anharmonic theory. Previous work has relied upon this approximation.liail3 

We may now clearly explain the anomalous diffuse scattering along the ridge between 
the Q = [222] bee Bragg peak and the N-point Q = |[552] position. This scattering is due 
to regions of the bee lattice fluctuating towards the hep lattice, with the result that intensity 
develops at the N-point, and shifts along the [112] direction. These are precisely the modes 
required to take the bcc lattice towards the hep lattice, as in Eq. |^. This interpretation is 
consistent with previous simulations of the bee lattice.By 

We have also examined the microstructure development during the transformation. In 
the beginning of the simulations, there are four equivalent (111) directions in the bee lattice. 
(Previous simulations have chosen periodic cells that break the symmetry between these 
directions. &ii) The final microstructure corresponds to domains of three different hep vari- 
ants, sharing a common (1120) direction along one of the original (111) directions. These 
domains are extended along this direction as well (corresponding to "needle-like" domains). 



7 



In addition, there is a domain in the metastable fee lattice. The final domain structure 
does not minimize the potential energy; a lower energy microstructure with two hep do- 
mains would satisfy the constraints imposed by the cubic simulation cell. Thus, the domain 
structure must arise out of the dynamics of the transformation process.0 

The hep variants are separated by {1011} twin boundaries, in accordance with exper- 
imental observations. These boundaries lie approximately along the bcc {112} planes, as 
shown in fig. ^ However, in some regions, these boundaries are not fiat, caused by basal 
plane stacking faults within the hep variants. Pinsook has discussed the formation of these 
stacking faults in terms of "transformation plasticity" ,0, and argues that these are neces- 
sary for the boundaries to lie parallel to the bee lattice's {112} planes, which intersect at 60° 
degrees, instead of the slightly different angles formed by an ideal {1011} hep twin boundary 
(55.9° for an ideal c/a ration). However, in our simulations, the boundaries often signifi- 
cantly deviate from being parallel to the bee {112} planes, indicating that more stacking 
faults occur than would be indicated by these geometric arguments. In both our simulations 
and those in ref. p^-p^, the hep basal plane stacking fault energy is lower than expected for 
Zr,ii which makes comparison with the experimental system difficult. 

The dynamics of the transformation correspond to a localized nucleation event, followed 
by growth of the domains. This is seen quite clearly in fig. ^ Such a picture is inconsistent 
with recent simulations of a dislocation-nucleated bcc—^hep transformation,!^ in which a bcc 
b = (100){011} dislocation nucleates the transition via a long wavelength distortion that 
evolves into a set of twin boundaries. However, this latter work is dependent on the particular 
dislocation; similar work with a b = l/2(110){li0} dislocation showed different dynamics.^ 
This work is also different from that of Pinsook and Ackland^, who observed plate-like 
geometry very quickly from a quenched simulation. In this latter work, the transition is 
homogeneous; unlike the work presented here, however, one of the (lll)^^^ directions is 
singled out in the simulation cell. The effect of this on the dynamics is unclear. They 
also observed an fee region, but in their simulations, this was a transient phase that was 
eliminated during the subsequent microstructural evolution. The presence of this region is 
likely to be sensitive to the choice of interatomic potential; for the work presented here (and, 
to a lesser extent, the potential used in Refs. |2^-plD, the difference between the cohesive 
energies of the fee and hep phases is small compared to the actual value in Ti and Zr, making 
the appearance here more likely. 

The present work provides a clear demonstration of the relationship between the dy- 
namics of the bcc^hep transformation, and the anomalous scattering observed in the par- 
ent phase. The observed zone-dependence does not require defects or static transformed 
(or distorted) regions in the lattice; this is clear from the simulations presented here and 
previously,^, and from experimental studies of pure metals. This is an important step in 
understanding the experimental results on the pre-transition dynamics, as measured through 
scattering experiments. Further work is required to understand how these dynamics are af- 
fected by defects, compositional fiuctuations in alloy systems, and the competing bee — > uj 
transformation. 
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FIG. 1. Potential energy versus time step, for the 18 x 18 x 18 and 36 x 36 x 36 system sizes. 
The large change in the potential energies is due to the bcc^hcp transition. 



FIG. 2. Final microstructure of 36 x 36 x 36 system, as viewed down the [111] direction, in a 
completely relaxed (zero force) configuration. The different atomic shadings indicate the potential 
energy of the atoms, with large, dark atoms having the highest energy; large, gray atoms having 
the lowest energy; and the smaller, lightest atoms having intermediate energies. This highlights the 
twin boundaries between hep variants. There are three such variants in the final microstructure, 
extending along the [111] direction, plus an fee region as indicated in the schematic picture shown 
in the lower right part of the figure. 



FIG. 3. Time-averaged positions of atoms during the steps indicated. The colors indicate the 
time-averaged potential energy, with blue being the lowest and red the highest. 

FIG. 4. The static structure facture 5(Q) calculated from simulations of the bee phase, as 
shown in the (110) scattering plane. The x-axis corresponds to the [110] direction; the y-axis 
corresponds to the [001] direction. 



FIG. 5. Time-dependent scattering at the 12 different Q = 3(774) wave vectors. As can be 
seen, those wave vectors whose directions are closest to [111] exhibit large scattering during the 
transformation. Note the differences in scales between graphs. 



FIG. 6. Time-dependent scattering at the 12 different Q = 3(558) wave vectors. Unlike 
the scattering at the Q = ^(774) wave vectors, there is little change in scattering during the 
transformation. 



FIG. 7. 'S'(Q) calculated from the time-averaged structures (as shown in fig. ^ for different 

time steps. The scattering plane is the same is in fig. ^. 

FIG. 8. <S'(Q) calculated from the time-averaged structures (as shown in fig. ^ for different 

time steps. The scattering plane is the same is in fig. |^. 



11 



This figure "figl.gif" is available in "gif" format from: 



http://arXiv.0rg/ps/cond-mat/0 1 03278 V 1 



This figure "fig2.gif" is available in "gif" format from: 



http://arXiv.0rg/ps/cond-mat/0 1 03278 V 1 



This figure "fig3.gif" is available in "gif" format from: 



http://arXiv.0rg/ps/cond-mat/0 1 03278 V 1 



This figure "fig4.gif" is available in "gif" format from: 



http://arXiv.0rg/ps/cond-mat/0 1 03278 V 1 



This figure "fig5.gif" is available in "gif" format from: 



http://arXiv.0rg/ps/cond-mat/0 1 03278 V 1 



This figure "fig6.gif" is available in "gif" format from: 



http://arXiv.0rg/ps/cond-mat/0 1 03278 V 1 



This figure "fig7.gif" is available in "gif" format from: 



http://arXiv.0rg/ps/cond-mat/0 1 03278 V 1 



This figure "fig8.gif" is available in "gif" format from: 



http://arXiv.0rg/ps/cond-mat/0 1 03278 V 1 



